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Abstract 

As the first useful Quantum Computers presumably will be quantum simulators, here 
the minimum number of qubits necessary for the solution of Schrodinger's equation 
in simple test problems is evaluated. From the present preliminary results it appears 
that it is possible to realize a useful quantum simulator with a register of only 10- 
15 qubits. The intrinsic sensitivity to some errors appears to be moderate without 
the need of error correcting methods. At present there is at least some indication 
that the amplitude errors are more dangerous than phase and decoherence errors. So 
corrections limited to amplitude errors may be very useful. 



1. Introduction 

After the pioneering work of Shor on the factorization [1] at present there is a 
growing interest in the study of quantum computers (QC) and the related problems 
of error control due to inaccuracies of the quantum gates, decoherence and so on (for a 
still useful review see [2]). As it is by now well known the QC may have an exponential 
advantage over the classical ones on many problems (for a very incomplete list see [3, 
4, 5, 6, 7]), but due to implementation difficulties the first QC will not have many 
qubits. So it is almost obvious that the first useful QCs will be quantum simulators 
and not true numerical computers. 

This is the first of a series of papers devoted to the simulation of a quantum 
computer for the solution of a generic quantum mechanical problem. The point of 
view is that of an experimentalist wanting to know how many qubits are necessary 
and what amount of errors it is possible to tolerate. So the present problem is to find 
the minimum requirements for the setting-up of a useful simulator to solve simple 
problems. 

The simulation of a quantum computer on a classical one was often done in the 
past [8], but for very specific problems. Here we study a system of general use to 
solve standard undergraduate problems in quantum mechanics as can be found in 
an exercise book [9]. Such a quantum simulator (QS) was simulated on a classical 
computer. Obviously the efficiency of such a simulation is exponentially slow, so 
nothing is gained in comparison with a direct solution of the various problems studied. 
From the experience gained so far, it is fair to say that a simulation of a QS on a 



classical computer approximately requires a computational power similar to the one 
necessary for the direct solution on the classical computer, except for the simulation 
of some errors. The simulation of the various kind of errors is the most uncertain 
part of this work, because at present there arc not sufficient experimental data on the 
sources of errors in the various implementations of real systems. It was then necessary 
to start with very simple hypotheses. Various problems were simulated: transmission 
and reflection of a Gaussian wavcpackct on various barriers, unidimensional harmonic 
and anharmonic oscillators, the solution of the radial equation of the tridimensional 
harmonic oscillator etc.. Here we report only the results about the unidimensional 
harmonic oscillator as the results of this simpler case are very similar to those of the 
other problems studied. 

2. Solution of the Schrodinger equation 

The algorithm we will use has been described elsewhere [10, 11] so let us merely 
show the essential points. The Schrodinger equation: 

ihj-^i^{x;t) = [-^^ + V{x)Mx;t) = [Ho + V{x)Mx;t) (1) 

may be easily solved on a QC provided that it is possible to decompose the algorithm 
in a sequence of unitary operations enough simple and, what is more important, 
efficiently implementables on a simple quantum computer. 

The first step is to decompose the region of interest of the space coordinate in 
2" intervals of length A and to represent these intervals in the Hilbert space of a n 
qubits register. This space is the direct product of n copies of the space of a single 
qubit: 

\in-i) <8 \in-2) O • • • <8) \io) = \i) (2) 

where (i„-i, in-2, io) is the binary decomposition of the index i. 

Then, neglecting the simbol (8), the wavefunction %p{x]t) is represented by the 
superposition: 



i^{x;t) = Ci{t)\i) = -E ^^"-.*"-.-.^oWKn-l>K„-2)-No) = 

i=0 Zn-i=0 in— 2=0 io=0 

2"-l 

= ^ E ^(^' = (^ + 0-5)A;t)N> (3) 

where A'^ is a normalization factor as the samples are not normalized. From this 
representation it is evident the exponential advantage to use a QS instead of a classical 
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computer: a register of n qubits can store 2" — 1 complex mimbers. This exponential 
advantage is completely lost in the simulation of the QS on a classical computer. 

To implement the various operations to solve the (1) it is necessary to extend the 
Hilbert space of the QS given by the register \n) with an ancella register \m)a of m 
qubits. So the complete Hilbert space of the QS now reads: 

\m)a <8) |n) (4) 

For the simulation of some errors of this QS on a classical computer it is neces- 
sary to add in the classical simulator the simulation of other qubits \l)i to describe 
unwanted quantum degrees of freedom or a rough representation of the environment. 
The complete Hilbert space of the classical simulation of the QS thus is: 

\l)i O \m)a O \n) (5) 



and for the classical simulator it is necessary to store 2('+™+") complex numbers. The 
experience gained shows that the practical limits of the simulation appears to be due 
more to problems of storage than to computational power. 

The Eq.(l) may be formally integrated by propagating the initial wavefunction 
every timestep e: 

ij{x;t + e) = e-^[-f^° + ^(^)l^ • (6) 
Using a time interval e small enough to neglect terms in it is possible to write: 

^-±[Ho + v{x)]e _ e"*^"" • e-*^(^)' (7) 

These operators are still unitary, simpler than the initial one and efficiently im- 

plementable on a QC. As the Fourier transform (JT) is also easily implemented [2]: 
by calling k the transformed variable, it is possible to write the first operator in the 
right hand side of (7) as: 



and so in the space representation: 

g-^ffoe = jr-i . (J . jr (9) 
The problem is now reduced to the implementation of an operator of the form: 



x) (10) 



This is possible using an ancella quantum register \m)a by means of the following 
steps: 

|0)a ® \X) ^ \f{x))a ® \X) ^ e'^(-) \f{x))a ® \x) ^ (11) 



3 



^ |0)a \x) = |0)a ® e^-^'(^) |x) (12) 

The first step is a standard binary representation of the function f{x) on the 
ancella register and may be implemented by means of c-not gates. The second step 
is the operation e""|m)a done on the ancella qnbits and so efficiently implemented 
by means of standard techniques with phase gates giving a phase shift to every qubit 
proportional to the weight of the qubit. The third step is just the reverse of the 
first one and immediately implemented by the same array of gates as the first one 
but applied in the reversed order. It may be useful to note that, in the examples 
studied, the implementation required hundreds of c-not gates. As the ancella qubits 
are retiirncd to the standard configuration it is possible to use the same ancella qubits 
for the kinetic energy and the potential energy computations. The potential energy 
part of the (7) is implemented in a similar way. In conclusion this algorithm can 
be implemented with no problems, except for some tricks necessary to minimize the 
number of the space and ancella qubits. 

3. Dimension of the registers 

As long as at present qubits are a very expensive resource, the first important 
problem to be solved is to determine the minimum number of qubits. Here we devote 
a special care to this problem as it is usually totally neglected. 

For the computation are needed n qubits to represent 2" points of the coordinate 
space and m ancellae qubits are also needed for the computation of the potential and 
the kinetic energy with the resolution of 2™ steps. The most simple determination 
is the number n of the 2" space points. This problem is common to every numerical 
representation of functions and is solved in a well known way by referring to the 
sampling theorem. As this theorem [12] is of central importance to keep the number 
of the qubits as small as possible it may be useful to remember some details. This 
theorem allows to compute the errors in reconstructing a function from its samples. It 
is enough to remember that the sampled function may be represented as the product 
of the original function times a train of Dirac delta functions. As long as the transform 
of a train of Dirac deltas is still a train of Dirac deltas, one gets replicas of the Fourier 
transform of the original function, shifted in frequency by the interval determined by 
the sampling frequency. So for a function whose spectrum is limited to a maximum 
frequency it is enough to sample at a frequency at least the double of this maximum 
frequency to avoid the superposition of the tails of the spectrum. At a lower sampling 
frequency there is an error in reconstruction due to tails superposition. This kind of 
error is referred to as aliasing error. To be concrete in the present problem it was 
assumed an initial wavefunction with a Gaussian shape. Since the Fourier transform 
is still a Gaussian so that the spectrum goes to zero very rapidly a surprisingly small 
number of samples was sufficient helping to keep the space register small. 

The implementation of functions is relatively simple: it is enough to "write" the 
values of the function in the binary representation in an ancella register of m qubits 
with a resolution of 2™ steps. This writing is realized by means of an array of c-not 
gates (the details will be given elsewhere) . A hard task is the evaluation of the number 
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m of the anccUa qubits necessary for the potential and kinetic energy computations. 
For this problem unfortunately there is not a theorem similar to the sampling theorem 
but only a less useful one on the characteristic function of the probability density of the 
function to be digitized (for a tutorial on these problems see [13] and cited references). 
So it was necessary to use some numerical tests. 

An unexpected result was that usually the number of the ancella qubits outweights 
the number of the qubits representing the space coordinate. The most sensitive part 
of this problem was the computation of the kinetic energy after the Fourier transform 
of the space wavefunction. 

A first trick to keep this number as small as possible is to tisc wavefunctions smooth 
in the space coordinates, in order not to have strong components of short wavelength. 
To have a concrete idea of what is happening, in problems with a simple potential as 
the harmonic oscillator in the case of an initial Gaussian wavepacket with a spread 
about 1/50 of the range of the space coordinate, and using only n = 6 qubits, i.e. 
2^ = 64 points for the space coordinate, the components of the Fourier transform 
of the wavefunction are populated only in the very lower part of the spectrum. In 
this case, for the computation of the kinetic energy up to to = 12 ancellae qubits 
were necessary to reach a sufficient resolution. A second trick to keep this number 
small was to compute by means of the ancellae the kinetic energy limited to the 
most populated part of the spectrum and to disregard the rest, hence using the full 
resolution of 2™ steps only for the used part of the spectrum. Tipical values used 
after various tests, were n = 6 space qubits and m = 7 ancella qubits. A total of only 
13 qubits was sufficient for most tests. 

4. Sensitivity to memory errors 

A further problem is the evaluation of the sensitivity of the simulator to errors. 

In the past various studies of the sensitivity to errors were done on numerical QCs 
[14, 15, 16, 17] and quantum computers dedicated to special problems [7], but here 
we concentrate on simulators able to solve a generic quantum mechanical problem. 
It is necessary to define exactly which kind of errors are to be studied and then to 
give an evaluation of the effect of these errors. The effect of the errors depends on 
the observable of interest. Here the effect was estimated by computing the fidelity 
defined by the comparison of the erroneous wavefunction |'!/')err with the correct one 
computed using the same system and zero errors. 

/ = M^)err\' (13) 

This quantity gives the probability that the erroneous wavefunction will pass a 
test for the correct one. The first kind of errors studied are the memory errors. These 
errors are simulated by applying a generic unitary matrix with random parameters 
to one qubit at a time of the register representing the space coordinate. This matrix: 

^ cos9 exp{i a} ■ sinO 

— exp{i P} • sinO exp{i {a + (3) ■ cos9} 

is applied at every timestep with random parameters with uniform distribution be- 



(14) 
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twcen plus and minus a given maximum (only a single parameter at a time is different 
from zero: as the errors are assumed to be small the probability of two or more si- 
multaneous errors is taken to be negligible). To ease the comparison of the effect of 
different errors the same string of random numbers was used in all the tests. To be 
more precise, the integration of Schrodinger's equation was done for 40 timesteps and 
the error was applied every step between steps 10 and 30. 

The dependence of the fidelity on the error parameter, the maximum amplitude of 
the error and the qubit in error have been explored. In the parametrization used the 
parameters a and /3 have an analogous effect so it is enough to study only the effect 
of the errors described by the a and 6 parameters. For brevity we refer as a, [3 and 6 
errors the errors due to the U matrix with only the specified parameter different from 
zero. Also for error ^ = 0.1 radiants we mean a series of U matrices with only the 
parameter 6 different from zero and random with uniform distribution bcetwen ±0.1 
radiants. Similarly for the other parameters. The results show a moderate effect of 
errors a (and so also and a strong effect of the 6 errors. 

The errors are applied from the steps 10 to 30 of the 40 integration, so that the 
fidelity starts with the value 1, remains in this value until the first application of 
the error matrix and then decreases randomly to the minimum when the last error is 
applied. In Table 1 we report this final value of the fidelity resulting from the error 
matrix applied to the various qubits of the register representing the space coordinate. 

Table 1. Final fidelity from a and 6 errors 



qubit 


a = 0.30 


e = 0.05 


e = 0.10 


e = 0.20 


e = 0.30 





0.90 


0.99 


0.96 


0.83 


0.66 


1 


0.90 


0.99 


0.96 


0.85 


0.70 


2 


0.90 


0.99 


0.97 


0.89 


0.76 


3 


0.83 


0.99 


0.96 


0.86 


0.71 


4 


0.96 


0.98 


0.91 


0.69 


0.43 


5 


0.98 


0.98 


0.93 


0.75 


0.53 



The first column indicates the qubit acted-on by the error matrix (14). The second 
column labeled (a = 0.30) gives the final fidelity obtained with errors of maximum 
0.30 radiants on the parameter a of the matrix (14). The following 4 columns give 
the final fidelity for errors with the parameter respectively of maximum 0.05, 0.10, 
0.20 and 0.30 radiants. By comparing the second and the last column it appears that 
the average effect of the errors in parameter a are almost negligible in comparison 
with the error in parameter 6. It is possible to note an unexpected low sensitivity to 
6 errors as high as 0.10 radiants, as shown in the 4** column. Another unexpected 
result is the low correlation of the fidelity with the number n of the qubit in error (of 
weight 2") indicated in the first column. 

The figures 1 and 2 give a plot of \xp{x]t)\'^ for errors in the qubits "0" and 
"5". It is useful to give a physical interpretation of these errors. This error gives a 
probability of reversal of the value of the qubit acted-on by the error. Upon indicating 
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with A the step in the space coordinate, the 9 error in the qubit "0" gives a reshuffling 
of the wavefunction values in adiacent intervals. A 6 error in the qubit " 1" gives an 
analogous reshuffling for intervals 2A, and so on for the 6 errors in the other qubits 
of higher weight. 

As the wavefunction used is a Gaussian wavepacket there are at least two different 
regimes depending on the width of the reshuffling interval, namely if it is smaller or 
greather than the width of the wavepacket. This effect is clearly evident in figure 1 
and 2, referring respectively to 9 errors in the qubits "0" and "5" and so referring to 
space intervals of lenght A and 32A respectively. 

5. Effect of the leak errors 

Other sources of errors are quantum leaks, gate errors and decohercncc. Here wc 
limit ourselves to leak errors. The quantum leaks are due to unwanted interactions 
of the qubits of the working register of the QS with quantum degrees of freedom not 
involved in the computation. To simulate these errors it is necessary to add some 
auxiliary qubits \l)i to the simulator Hilbert space. So this space now is given by (5). 

Here we study the effect of the leak errors due to the interaction of a single qubit 
of the space register with a single extra quantum degree of freedom described by a 
single qubit starting in a pure state. So it is possible to represent this interaction 
by means of a unitary 4x4 matrix applied to a two qubit space given by the qubit 
supposed to be acted by the error and the leak qubit at the specified times and tracing 
away the ancella leak qubit at the end of the computation. Also in this case the error 
matrix is applied every timestep between steps 10 and 30. The error matrix is now 
4x4 and so it is in general described by too many parameters to be useful. So in the 
absence of sufficient experimental data we did select two particular very interesting 
matrices U\ and U2 depending from a single parameter given by [18]: 



Ux = 



10 

cos9 — sin0 

10 

sin6 cos6 



U2 



10 

cos9 sin9 

sin9 — cos9 

10 



(15) 



To have an idea of the reason of the choice of these matrices, it is useful to note 
that the matrix Ui, when applied to a system composed of a single qubit and a single 
ancella qubit gives a pure decoherence. In fact denoting with X, Y, Z the parameters 
of the Poincare sphere of the density matrix of a single qubit: 



1 + Z 
X + iY 



X - iV 
1 - Z 



(16) 



the matrix Ui applied to a system composed of a single qubit with this density matrix 
and an ancella qubit in the pure state |0) after the trace of the ancella qubit gives for 
the system reduced density matrix the transformation: 



X' = cos9X 



y = cos9Y 



Z' = Z 



(17) 



7 



and so the coherence {X — iY) is reduced. The analogous computation with the 
matrix U2 gives: 

X' = cose X Y' = COS0 Y Z' = sin^e + cos^O Z (18) 

and so the coherence is reduced and there is also a shift and a factor lower than 1 
of the Z component of the reduced density matrix representing the population. The 
fidelity resulting from errors due to the application of the Ui and the U2 matrices 
to the various qubits of the space register was computed in a way similar to the 
proceeding case but with the reduced matrix obtained by tracing away the ancella 
qubit 1^);. 

Similar to the previous case, the fidelity starts with the value 1 and then randomly 

decreases to the final vahie when the last error matrix is applied. The parameter 9 
is supposed to be random with uniform distribution between plus and minus the 
maximum value indicated and the error matrices are applied from the timestep 10 to 
the timestep 30. 



Table 2. Final fidelity from leak errors Ui 
quUt 6* = 0.05 61 = 0.10 6» = 0.30 






1.0 


0.99 


0.82 


1 


1.0 


0.98 


0.89 


2 


1.0 


0.99 


0.79 


3 


1.0 


0.99 


0.77 


4 


1.0 


1.0 


0.95 


5 


1.0 


1.0 


0.93 



Table 3. Final fidelity from leak errors U2 



qubit 


6 = 0.05 


9 = 0.10 


9 = 0.30 





1.0 


0.98 


0.76 


1 


1.0 


0.99 


0.83 


2 


1.0 


0.97 


0.82 


3 


1.0 


0.94 


0.78 


4 


0.99 


0.97 


0.79 


5 


1.0 


0.97 


0.66 



In Table 2 we report the final fidelity for leak errors on the various qubits of the 
register of space coordinate. The first column indicates the qubit acted-on by the leak. 
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The following three columns give the final fidelity for various values of the maximum 
of the error parameter 9 in the matrix Ui. The Table 3 give the analogous quantity 
for the matrix U2- In both cases there is a low correlation of the fidelity with the 
weight of the qubit in error. 

The leak error is due to the transfer of some of the quantum information on 
the system to a further quantum degree of freedom. But as the computation of 
the fidelity is done by neglecting this degree of freedom and calculating the reduced 
density matrix, part of the quantum information is regained. This effect may explain 
the (lower than expected) effect of the quantum leak, even with an error as high as 
the 0.30 radiants of the parameter used in the matrices Ui and U2- 

In figure 3 and 4 we report a plot of \tp{x;t)\'^ for the errors due to the matrix Ui 
applied to the space qubit "0" and "5" for the error ^ = 0.1 radiants, corresponding 
to the third column of Table 2. It is possible to note a very small effect of this kind 
of error. In figure 5 and 6 are reported the analogous results for the matrix U2- 
The matrix U2 generates, along with decoherence errors, also population errors. In 
this case it is possible to distinguish at least two different regimes. This last kind of 
error gives an effect similar to the error from the matrix U shown in Table 1 and 
figure 1 and 2. As in the case of the figure 1, in figure 5 it is evident a granularity 
due to the reshuffling of the wavefunction values in adjacent space intervals. This 
granularity changes to a random noise in figure 6 referring to errors in qubit "5" as 
in the corresponding U case reported in figure 2. 

By comparing the results obtained from the matrix Ui (decoherence) and those 
from the matrix U2 (decoherence and population), it appears that the last kind of 
errors is more dangerous. 

6. Conclusions and future developments 

The first aim of this work was to have reliable data on how many qubits arc neces- 
sary to implement a QS for the solution of a generic and simple quantum mechanical 
problem. Susprisingly the number of ancella qubits necessary for the implementation 
of the various operators outweights the number of the qubits representing the space 
degree of freedom of the system. But for the simple problems studied here, with some 
care a total of only 13 qubits proved to be sufficient. Another useful result is that the 
QS appears to be quite robust against the explored errors: an error of 0.1 radiants 
on all the parameters indicated gives a fidelity higher than about 0.9. Obviously for 
other observables and for a different measure of the error this conclusion may be not 
valid. The fidelity is largely independent from the weight of the qubit in error. 

To ease the comparison of the memory errors as defined here with other definitions, 
it is sufficient to say that the a and f3 errors are similar to the phase errors and the 
9 error is similar to the amplitude error. So it is possible to say that in this QS the 
amplitude errors are more dangerous than the phase errors. 

The most important development is to enlarge the single leak qubit to a larger 
space in order to model the environment. The interaction of a single qubit with the 
environment has been modeled by various authors[19, 20]; it is known that two ancella 
qubits are suHicient, but one needs to trace away these qubits at every integration 
step in order to wash-out any memory of the system taken by the environment. Hence 
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it is necessary to describe the system by means of the density matrix. This does not 
create any basic problem, but only problems of storage and computational power of 
the classical computer. 

Acknowledgments 

This work was supported by M.P.I. I would like to thank prof. V.G. Benza for his 
careful reading of the manuscript and the many useful suggestions. 



(*) E-mail giuliano.strini@mi.infn.it 
References 

[1] P. W. Shor. Proc. 35th Ann. Symp. Found. Comp. Sci. IEEE Comp. Soc. 
Pr. 124 (1994). 

[2] A. Eckert and R. Jozsa. Reviews of Modern Physics 68, 733 (1996). 
[3] S. Lloyd. Science 273, 1073 (1996). 

[4] D. S. Abrams and S. Lloyd. Phys. Rev. Lett. 79, 2586 (1997). 
[5] D. A. Lidar and H. Wang. Phys. Rev. E59 2429 (1999). 
[6] D. A. Lidar and O. Biham. Phys. Rev. E56 3661 (1997). 



[7] G. Benenti, G. Casati, S. Montangero and D. L. Shepelyansky. [arXiv: quant - 
ph/0107036| . 



J. Wallace Quantum Computer Simulators. A Review, http: / / www. dcs. ex. ac.uk Wallace 



[9] I. I. Gol'dman et al.. Problems in Quantum Mechanics. Infosearch Limited. 
London 1960. 

[10] C. Zalka. Fortschr. Phys. 46, 877 (1988). 
[11] S. Wiesner. ^rXiv:quant-ph/9603028 . 



[12] F. M. Reza. An introduction to information theory. McGraw-Hill 1961. 

[13] A. Grassi and G. Strini. Alta frequenza XXXIII E547 (1964). 

[14] K. M. Obenland and A. M. Despain. Proc. of High-Performance Computing 
HPC'98 228 (1998). 



10 



[15] C. Miquel, J. P. Paz and R. Perazzo. Phys. Rev. A54 2605 (1996). 

[16] A. Barenco and A. Ekert. Phys. Rev. A54 139 (1996). 

[17] C. Miquel, J. P. Paz and W. H. Zurek. Phys. Rev. Lett. 78 3971 (1997). 

[18] Chi-Sheng Niu and R. Griffiths. Phys. Rev. A60, 2764 (1999). 

[19] B. Schumacher. Phys. Rev. A54, 2614 (1966). 

[20] G. Strini. Unpublished Lecture Notes (2000). 



11 



Figure captions 



Figure 1. Plot of \'tjj(x;t)\'^ obtained with the memory error matrix U applied to 
the qubit "0" of weight A. The granulosity due to the reshuffling of the wavefunction 
values at a distance A is clearly visible. 

Figure 2. Plot of \'4'{x;t)\^ obtained with the memory error matrix U applied to 
the qubit "5" of weight 32 A. The reshuffling of the wavefunction values at a distance 
32A originates a random noise. 

Figure 3. Plot of \tp{x;t)\'^ obtained with the leak error matrix Ui applied to the 
qubit "0" of weight A. This leak originates a decoherence error and a very low noise. 

Figure 4. Plot of \tp{x;t)\^ obtained with the leak error matrix Ui applied to 
the qubit "5" of weight 32 A. As in the preceeding figure 3 this leak originates a 
decoherence error and a very low noise. 

Figure 5. Plot of \^{x;t)\'^ obtained with the leak error matrix U2 applied to the 
qubit "0" of weight A. This leak originates a decoherence error and an amplitude 
error. The amplitude error in the qubit "0" originates a reshuffling of the wavefunction 
values at a distance A and so a granulosity very similar to the one due to memory 
errors reported in figure 1. 

Figure 6. Plot of \il){x;t)\^ obtained with the leak error matrix U2 applied to the 
qubit "5" of weight 32 A. This leak originates a decoherence error and an amplitude 
error. The amplitude error in the qubit "5" originates a reshuffling of the wavefunction 
values at a distance 32A and so originates a random noise similar to the one due to 
memory errors reported in figure 2. 
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